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Abstract 

The equation of state of a system of fermions in a uniform magnetic field is obtained in terms 
of the thermodynamic quantities of the theory by using functional methods. It is shown that the 
breaking of the 0(3) rotational symmetry by the magnetic field results in a pressure anisotropy, 
which leads to the distinction between longitudinal- and transverse-to-the-field pressures. A crite- 
rion to find the threshold field at which the asymmetric regime becomes significant is discussed. 
This threshold magnetic field is shown to be the same as the one required for the pure field con- 
tribution to the energy and pressures to be of the same order as the matter contribution. A 
graphical representation of the field-dependent anisotropic equation of state of the fermion system 
is given. Estimates of the upper limit for the inner magnetic field in self-bound stars, as well as 
in gravitationally bound stars with inhomogeneous distributions of mass and magnetic fields, are 
also found. 
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I. INTRODUCTION 



The largest magnetic fields observed in nature are associated to some of the most extreme 
astrophysical objects, the compact stars. For pulsars the typical magnitudes of surface 
magnetic fields have range ~ 10^^ — 10^'^ G [1]. The measured periods and spin down of soft- 
gamma repeaters (SGR) and anomalous X-ray pulsars (AXP), as well as the observed X-ray 
luminosities of AXP, indicate that a certain class of neutron stars called magnetars can have 
even larger magnetic fields, reaching surface values as large as 10^^ — 10^^ G [2]. Furthermore, 
if the suggestion 3|] that these stars can be the source of gamma-ray bursts is confirmed, 
their magnetic fields should be even larger in order to drive a substantial Poynting fiux- 
dominated relativistic outflow. Up to now, about ten highly magnetized neutron stars have 
been identified in our galaxy, but based on the population statistics of SGR it is expected 
that magnetars constitute about 10 % of the neutron star population 

The existence of the strongest magnetic fields in compact stars poses the question of their 
origin. The simple hypothesis that a relative small magnetic field of a progenitor star can 
be amplified during the star's gravitational collapse due to magnetic flux conservation j^, 
cannot even substantiate the high values of the surface fields in magnetars Gj. Another 
generation mechanism is the so called magnetohydrodynamic dynamo mechanism (MDM). 
The MDM is based on the amplification of a seed magnetic field due to the rapidly rotating 
plasma of a protoneutron star. It is generally accepted as the standard explanation for the 
origin of the magnetar's large magnetic fields. For the MDM to explain the large surface 
field strengths observed in magnetars, the rotational period of the protoneutron stars that 
originate them should be < 3ms. Nevertheless, this mechanism cannot substantiate all the 
featu..es of the supernova .e,„„a,.s ...ounding these ob.eets B. Pa. of the rotational 
energy is supposed to power the supernova through rapid magnetic braking, from where it 
is inferred that the supernova remnants associated with magnetars should be an order of 
magnitude more energetic than typical supernova remnants. However, recent calculations 
[3] indicate that their energies are similar. In addition, one would expect that when a 
magnetar spins down, the rotational energy output should go into a magnetized particle 
wind of ultrarelativistic electrons and positrons that radiate via synchrotron emission across 
the electromagnetic spectrum. Yet, so far nobody has detected the expected luminous pulsar 
wind nebulae around magnetars ^ . On the other hand, some magnetars emit repeated flares 



or bursts of energy in the range of 10^^ — W^^erg 10|]. Since the emitted energy significantly 
exceeds the rotational energy loss in the same period, it is natural to expect that the energy 
unbalance could be supplied by the stellar magnetic field, which is the only known additional 
energy source. Nonetheless, from the spin history of these objects, there is no clear evidence 
of any surface magnetic field damping [1 Ij . 

From the previous considerations it is clear that alternative mechanisms to the standard 
magnetar model [2i] should be explored. A reasonable approach is to investigate possible 
microscopic mechanisms, based on the quantum phase of the core, through which a seed 
inner star's magnetic field can be generated and/or boosted. Some propositions in this 
direction already exist in the literature jo], [l2|. Even though, to find a connection between 
the microscopic phase of the star's core and the astrophysical observations, other important 
properties of the star's core matter need first to be better understood. Along these lines, a 
very important problem to elucidate is the infiuence of a magnetic field on the star's equation 
of state (EoS). 

Over the years, many works have been dedicated to the effects of magnetic fields in 

. However, in general, 

when finding the field-dependent contributions to the energy density and pressures, they did 
not follow a unique and consistent scheme, thereby different papers have different stands on 
what should be the correct field contributions to the pressure and energy. Moreover, it is 
known that in the presence of a magnetic field the pressure splits in two terms: transverse 
and parallel to the field direction, due to the breaking of the spational rotational symmetry. 
Nevertheless, some authors ignored the pressure anisotropy even at very strong magnetic 
fields, where it becomes significant. Besides, one can identify, depending on their origin, two 
different field contributions to the energy and pressures. One coming from the magnetized 
matter and the other from the Maxwell term. Despite all this, some of the previous studies 
do not take into account the pure field effect (coming from the Maxwell term) in the energy 
density and pressures, even though it is always present and in some limits it can be the 
dominant one. 

The main purpose of this paper is to develop a systematic and self-consistent approach to 
treat the EoS of a magnetized system. We analyze under what conditions the pure magnetic 
contribution to the energy and pressures is much smaller than the matter contribution, 
as well as when it is self-consistent to neglect the differences between the transverse and 



parallel pressures (isotropic limit). We carry out our study in a theory of free fermions 
only interacting with an applied uniform and constant magnetic field, but the method we 
developed to analyze the effect of the different contributions to the EoS can easily and 
straightforwardly incorporate interactions. 

Another outcome of our investigation is an improved estimate for the upper limit of com- 
pact stars' inner magnetic fields. Previous estimates were done assuming a gravitationally 
bound star with spherical and homogeneous mass distribution and a uniform magnetic field. 
For gravitationally bound stars we demonstrate that when the homogeneous mass density 
and uniform field distribution conditions are relaxed, the inner field in the high-dense core 
can attain values two orders of magnitude larger than previously found. We also estimate 
the inner magnetic fields in self-bound stars, which results of the same order than that of 
inhomogeneous gravitational bound stars. Using our magnetized fermion model, we calcu- 
late the threshold field that separates the isotropic and anisotropic regimes. This field turns 
out to be smaller than the estimated upper values of the stars' inner fields, indicating that 
the anisotropic effects can be relevant for the physics of the cores. 

The outhne of the paper is as follows. In Sec. II we estimate the upper hmits for the 
magnetic field in self-bound and gravitationally bound compact stars. The derivation of the 
Maxwell and Dirac field contributions to the covariant energy-momentum tensor is reviewed 
in Sec. III. In Sec. IV, the quantum-statistical average of the energy-momentum tensor 
components are calculated using a functional method. From these results the system energy 
density, and the parallel and transverse pressures are obtained in terms of the thermody- 
namical quantities. A covariant structure for the energy-momentum tensor, in agreement 
with the symmetries of the magnetized many-particle system, is given in terms of the ther- 
modynamic quantities. In Sec. V, the EoS of the magnetized fermion system is found at 
zero temperature and finite density. Numerical results for the energy density and pressures 
as functions of the magnetic field are presented and the significance of the matter and field 
contributions for the different ranges of densities and magnetic field strengths are discussed. 
Also, the threshold field for the transition between the isotropic and anisotropic pressure 
regimes is obtained. Our concluding remarks are stated in Sec. VI. In the Appendix the 
system thermodynamic potential is derived using Ritus's eigenfunction method. 
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II. COMPACT STAR'S FIELD ESTIMATES 



Since the interior magnetic fields of neutron stars are not directly accessible to observa- 
tion, their possible values can only be estimated with heuristic methods. A direct application 
of the virial theorem leads to inner field estimates of order ~ 10^® G [3] for compact stars 
with masses M ~ 1.4Mq and radius R ~ 10~^Rq. This derivation was done for gravita- 
tionally bound objects with uniform fields and mass density. In this Section, we show that 
if these conditions are relaxed the star's inner magnetic field may reach even higher values. 
Below we consider two possibilities: 1) self-bound objects with uniform magnetic fields, and 
2) gravitationally bound objects with a physically more realistic case of inhomogeneous field 
and mass distributions. 



A. Self-bound compact stars 

Self-bound stars are stars made of stable u, d and s quark matter. Let us explain how 
this can be possible at least from a theoretical point of view. If the star's density is high 
enough for deconfinement, a quark matter phase may occur. Under these conditions, the 
up and down quarks can convert into s-quarks via weak interactions. In fact, the quark 
system will prefer to do so in order to lower the Fermi energy by increasing the degeneracy. 
The charm, top and bottom quarks are not relevant in this analysis because their masses 
are much larger than the strange's and the typical densities of the stars will not be enough 
to produce these other flavors. The thereby formed three-flavor quark matter, composed 
of a mix of m, d and s quarks, is known as strange quark matter. It has been conjectured 
{igI . jl?! that at zero pressure the strange quark matter will have a lower energy per baryon 
than ordinary matter, which has e Matter = 939 MeV/baryon. This possibility would make 
the strange matter the most stable substance in nature. According to this, ordinary nuclei 
would lower their energy by converting to strange matter, but it has been shown that the 
conversion rate is negligible under almost all conditions, except perhaps in neutron stars 

n nn 

18[. Therefore, the possible existence of strange stars cannot be ruled out [17|, |18| . 

Assuming there are strange stars out there, a reasonable question to ask is: how big 
a magnetic field can be sustained by them? Energy-conservation arguments can help to 
estimate the maximum field strength; one expects that the magnetic energy density should 



not exceed the energy density of the self-bound quark matter. Based on this, the maximum 
field allowed is estimated as 

B_ . < (52^ . 1^5 X 10-G. (1) 

enc ehc 

We call attention that it results two orders of magnitude larger than the estimates done for 
gravitationally bound stars assuming uniform field and mass density. 



B. Gravitationally bound stars with inhomogeneous mass and field distributions 

Let us consider now the case of gravitationally bound stars. As known, neutron stars 
are the remnant of type-II supernova explosions. In a neutron star, pressure rises from zero 
(at the surface) to an unknown large value in the center. Also the density changes from 
surface values much smaller than the saturation density, ~ 4 x 10^^ g/cm? -the density 
at which nuclei begin to touch- to inner values several times the normal nuclear density 
Pn = 2.8 X 10^^ gjcvn?. At such high densities, deconfinement can occur and the star's 
core can have quark matter. Thus, the cores of very massive neutron stars are the best 
natural candidates for the realization of the transition from hadronic matter to a deconfined 
quark phase. This possibility was pointed out by several authors long time ago [l9^. Stars 
with radius- dependent density leading to confined nuclear matter in the outer region and 
deconfined quark matter in the core are called hybrid stars. For excellent reviews on this 



topic see [20 1. 



What can we say about the star's magnetic field in this case? Well, we know that the 
stellar medium has a very high electrical conductivity, hence the magnetic flux is conserved. 
Since the flux is conserved, during the formation of the neutron star, that is, during the 
protoneutron star period, the magnetic field strength should increase with the increase of 
the matter density. In addition, some of the proposed phases [3l-j2l| that could be realized 
in the core of the compact object that result at the end of the supernova explosion might 
also contribute to increase the field strength even more. For instance, if the resulting core is 
dense enough as to be in the color superconducting color-flavor-locked (CFL) phase and a 
magnetic field of strength comparable to the Meissner mass of the "rotated" charged gluons 



were present, this field would trigger an instability 2l| which in turn would lead to the 



generation of a gluon vortex condensate. The gluon vortex state so produced behaves as 



a paramagnet, thus increasing the total magnetic field in the core. This vortex phase is 



known in the literature as the paramagnetic CFL (PMCFL) phase 2l|. As these arguments 



indicate, the magnetic field in the core could be higher than in the surface. 



Let us briefiy outline the field estimate obtained in Ref. 



15| for the case of homogenous 



field and mass distributions. From the equipartition theorem under these conditions 

I^UA ^ = '-G^, (2) 

3 / Svr 5 i? ' ^ ^ 

they found that the maximum field was given by 

^"^'^^ = ^® (^) (:^) ' 

Using Hq = 2 X 10® G, and taking into account that a typical neutron star has M ~ IAMq 
and R ~ (0.14 x 10~'^)Rq, one easily estimates the maximum strength as Hmax ~ 10^® G. 

But as we said above, we are interested in the more realistic situation where the con- 
straints of uniform mass and field distributions are relaxed. With this aim in mind, let us 
assume that both the mass density and the magnetic field increase from the surface (r = R) 
to the star center (r = 0) and let us consider the following parametrizations for the density 
and magnetic field respectively 



P = Po 



a > 0, (4) 



H{r) = Hs + {Ho - Hs) b > 0, (5) 

Po is the density at the core; and Hq are the magnetic fields at the surface and inner core 
respectively; and a and b are parameters to be determined. 

Using the mass density (HI), a star with a spherical configuration of radius R will have 
mass 



The coefficient a must be positive, but apart from that it is totally arbitrary. We can use 
typical values of neutron stars' mass M = 1.4M0 and radius R = 10km, as well as realistic 
core density estimates to find a region of physically acceptable values for a. The results 
are shown in Fig.l-(a). The plot gives the mass density coefficient a as a function of the 
parameter n that characterizes how much larger than the nuclear density is the star's core 
density po = np^, n = 1,2, .... Notice that n must be larger than 2 to obtain an acceptable 




FIG. 1. (a) The mass density coefficient a vs the density parameter n. (b) The core magnetic field 
Hq vs the field coefficient b. 



(positive) value of a. For n = 3, a realistic core density case, the parameter a is positive 
and lies between 10 and 20. 

As is a, the parameter h is also arbitrary. To obtain a reasonable set of values for h we 
need to use the equipartition theorem in our inhomogeneous star model. The gravitational 
energy of a spherical distribution of mass with density (jl]) and radius R is 



Wg = 2n f p{r)(f){rydr = -An'^GplR^Fia) 

J 



with F{a) given by 



Fia) 



8a^ + 60a2 + STa^ 



15(a + 2)(a + 3)(a + 5)(2a + 5) 
The gravitational potential (j){r),r < i? in (J?!) is 

p{r')dV' 



I /j^ ~~rj^^ I 



1 / r \ ' 



: + 



a+2 



a + 2 3 Vi?/ (a + 2)(a + 3) Vi?. 
On the other hand, the magnetic energy corresponding to the field configuration 



R HHr] 



r dr = - 



1 



in 



2 Jo 



r'^dr 



3(6 + 3)(26 + 3) " 6(6 + 3)(26 + 3) 



tt2 

2(26 + 3) \ 



(7) 



(8) 



(9) 



IS 



(10) 
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Taking into account that magnetars' surface fields are ^ Hg = lO^^G, and considering 
the value of a that corresponds to = 3f)Af 0, we can use the equipartition of the magnetic 
(ITU]) and gravitational ([7]) energies, Wg = Wm, to graphically find the inner field Hq as a 
function of the parameter b. The resulting curve is plotted in Fig. l-(b). The smaller b is, 
the slower the field decays from the inner region to the surface, and thus the larger the inner 
field required to produce Hg = lO^'^G in the surface. For a value of b between 0.1 and 0.01 
the core field can be larger than (10^^ G). 

The above estimate for the field is model-dependent. Even though we make no claim that 
this ad-hoc model correctly describes the real way the field varies with the radius in a hybrid 
star, our derivation serves to illustrate how the more realistic assumption of inhomogeneity 
can be consistent with stronger field estimates in the core than those previously found in 
the literature using homogeneous distributions. Since at present there is no reliable way to 
know the exact dependence of the density and field with the radius in a real neutron star, 
we will have to wait for more observations to validate or not this possibility. 

We will see now, based on energy conservation arguments within a microscopic analysis, 
that there is a natural scale for the core magnetic field ~ 10^° G. Given the energy density 
per baryon at the core of the gravitationally bound system 

e P\\ N 

— = -— + , 11 

ua ua riA 

where the baryon number ua = + + Us) = N/3, let us assume that there will be 
some field value from where the parallel pressure becomes negligible {p\\ ~ 0) and then let us 
estimate the order of this maximum field just by reasoning that the magnetic energy density 
should be at most as large as the energy density of the baryon system. Then, neglecting the 
first term in the RHS of Eq. flTTl) . one has 

— = 3/i (12) 

Ua 

and consequently, 

H ~ 9/iVe^c (13) 

For the phenomenologically acceptable value of /i ^ 400 MeV, we get from f|T3|) that the 
magnetic field ~ lO^^G. As will be shown later in the paper, the parallel pressure indeed 



^ Densities of this order are large enough for deconfinement. 
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decreases with the field and for a system of free fermions becomes neghgibly small at a field 
strength of order lO^^G. 

Thus, we conclude that field strength of 10^'^ G makes a natural maximum scale for both 
self-bound and at the core of high-density inhomogeneous gravitational-bound compact stars. 



III. ENERGY-MOMENTUM TENSOR AT H^^ 

For the sake of completeness and understanding, let us outline the steps that lead to a 
symmetric and gauge invariant energy-momentum tensor of matter and fields. The energy- 
momentum tensor T^j^ is associated with the Noether's currents of an infinitesimal space- 
time-dependent translation and in terms of the Lagrangian density c^^V^) of the theory 
it is given by 

dC 

This T'^^ tensor is the so-called canonical energy-momentum tensor. In general it is neither 
symmetric non gauge invariant. Nevertheless, these two shortcomings can be fixed by taking 
advantage of the arbitrariness of T^'^, which is defined up to a total derivative 

T^""" T^"" + dpM^"" (15) 

with M^'^'^ = — M'"^^. The redefinition (fT5|) preserves the conservation {d^T^^ = 0) and the 
value of the global energy-momentum four- vector 

P^ = J d^xT^"^ = J (fxiTf"^ + dpM^^P), (16) 

as long as M**^^ vanishes sufficiently rapidly at infinity. Then, the energy and momentum, 
which are measurable quantities, as well as their conservations, are maintained despite the 
addition of a new tensor to the canonical one. With a suitable choice of the tensor M'^"'^ 
the energy-momentum tensor can be converted into a symmetric and gauge invariant one. 

On the other hand, there is a general derivation which guarantees from the beginning the 
symmetry and gauge invariance of T'^" (denoted by t^'^ from now on). The idea derives from 
the fact that if the matter fields are coupled to gravity, the energy-momentum tensor plays 
the role of the source of the gravitational field. In this case, the introduction of gravity will 
generalize the space-time transformations to the frame of the general covariance of general 
relativity, rather than to the particular Lorentz transformations. Thus, the program to follow 
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is to consider the starting theory in a curved space-time geometry where the Minkowski 
metric r^^i/ is replaced with the general metric g^^, the volume element d^x by d^x^—g 
(with g being the determinant of the metric), and then to obtain the energy- momentum 
tensor through the invariance of the action with respect to the variation of the metric 22 1 



i / d'x^dg.^T^"' = (17) 



2 

Once a manifestly symmetric and gauge-invariant t^^ tensor is obtained by this procedure 
we can switch off the gravitational field by returning to the Minkowski metric {g^^ — )■ r^^i^). 

To find the energy-momentum tensor of a fermion system in the presence of a magnetic 
field, we start from the Lagrangian density 

£(7/>, F^,) = £^,(F^,) + £^(^, F,,) (18) 

where £^^(F^y) denotes the Maxwell Lagrangian density and C^^ip, F^,^) that for the Dirac 
field in the presence of the external magnetic field. They are respectively given by 

Cj,^{F,,) = --F.^F'^^ (19) 



and 



£^(^, F,,) = - m)^ + ^H'd.Y - m)^ (20) 



with the right and left gauge covariant derivatives given respectively by 

= it,, - eA^ (21) 

and 

= -^^M - (22) 
where is the electromagnetic potential associated with the external applied field. 

A. Energy-momentum tensor of the Maxwell field 

To find the manifestly symmetric and gauge invariant energy-momentum tensor for the 
Maxwell field we have from ( IT7|) that 

- — ^(v^^aJ (23) 



-9 Sg, 



pX 
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where La^ is obtained from (fT^ by explicitly introducing the dependence on the metric 
tensor 

= —F.^F^.g^^g'^^ (24) 



Taking into account that 



5^ 1 



-99'' (25) 



S9px 



-9''9"' (26) 



one finds 



= -F^pp; - Ca^" (27) 

Returning to the Minkowski space and considering in particular the case for a constant 
and uniform magnetic field H along the xs-direction we have from ( 1271) that 

T^^ = {e- p)u^u'' + p(r7f - r/f ) (28) 



where e = p = iy^/2, is the medium four- velocity which in the rest system takes the 
value = (l,"(f), r^ii"^ is the longitudinal Minkowskian metric tensor with /i, i/ = 0,3 
and 7]^ is the transverse Minkowskian metric tensor with /i, z/ = 1,2. The fact that the 
energy-momentum tensor of the magnetized space becomes anisotropic, having different 
pressures in the longitudinal p\\ and transverse p_L directions {—p\\ = p± = is due 

to the breaking of the rotational symmetry 0(3) produced by the external field. As a 
consequence, the Minkowskian metric splits in two structures, one transverse rj'^ = F'^'^F^ 
(where F'^'' = F^^/ H denotes the normalized electromagnetic strength tensor) and another 
longitudinal r]f = r?^'' - F^'PF^. 



B. Energy-momentum tensor of the Dirac field 

To follow the previous approach in the case of spinor fields is more involved. The problem 
is that there is no representation of the GL{A) group of general relativity which behaves 
like spinors under the Lorentz subgroup [23|. Then, to put fermions in interaction with 
a gravitational field a new formalism is required. The first formulation of spinor fields in 
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:liemannian space-time was done in [2J] by introducing the so called vierbein or tetrad fields 
25|. The vierbeins V^{x) connect the Minkoskian metric r/^i, with the metric tensor of a 
general coordinate system g^p by the map 

9ap = V,^ix)V/^ix)v,. (29) 

where 



Kr(x) = ( ^ ) (30) 



x=X 



with Ux being the normal coordinates of a local Minkowski space at point X, and x" the 
corresponding general coordinates at that point. In this way, the geometries in general 
relativity can be described in terms of a vierbein field instead of the usual metric tensor field. 
We have adopted the convection that the indices given by Greek letters at the beginning 
of the alphabet (i.e. a, (3, etc.), are related to magnitudes in the general reference system, 
while those at the end (i.e. /i, z/, p, etc.) to magnitudes in the local Minkowsky reference 
system. 

On the other hand, the introduction of the vierbein fields also make it possible to gener- 
alize the algebra of the Dirac gamma- matrices 7^^, given by the anticommutation relation 

{r,Y} = 2v'''I (31) 

to the curved space 



{7^,7^} = 2^7"^^/ (32) 



with 7" = 1/"(x)7^. 

The Lagrangian density of the Dirac fields f l20|) can be rewritten in curved space-time 
with the aid of the vierbeins fields as 



1 r- 



\pYV^V^V,ilj - (V^tfj) V^V^Y^] - m^iP (33) 



where Vjy = VJ* (Da + r„) is the covariant derivative in curved space with connection Ta{x) = 
^T.'^'^V^ {x){-^Vf3,y{x)) and T.^'^ = ^ [7'^,7'^] being he generator of the Lorentz group. 
Taking into account that 

= det[V{x)], 5g^p = - {g^,V^ + gp,V^) (34) 
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it is found that the energy- momentum tensor obtained from (l33l) is 26 1 

V'^ix) 6det[V{x)]C^ 



det[V{x)] 5V^ 



(35) 



Returning to Minkowski space (i.e. with the replacement 7" — )■ •y^, — )■ -D^) we obtain 



(36) 



Let's assume a uniform and constant magnetic field H along the Xs-direction and use the 



Landau gauge A^^^ = Hxi5^j,2 in the covariant derivative = + ieA^j^. Then, we find 



1 
2 



(37) 



4%x) = tpzd^y'ij + C^, (38) 

and 

4\x)=t'^{D^y)^ + C^, J = 1,2 (39) 

Once again, the asymmetry between the longitudinal and transverse diagonal components 
of r^*^, is related to the breaking of the 0(3) symmetry. Moreover, since the Landau gauge 
is not symmetric, we have that although a magnetic field along the xs-direction conserves 
the rotation group 0(2) in the corresponding perpendicular plane, the asymmetry of the 
potential introduces an apparent asymmetry in the transverse indices (i.e. r^^ 7^ '^^^)- 
This apparent asymmetry, as we will see in the following section, cannot be present in the 
quantum-statistical average of r^*^, which can only depends on the field and not on the 
potential. 



IV. ENERGY AND PRESSURES OF THE DENSE MAGNETIZED SYSTEM 

As is known |i22|], in the reference frame comoving with the many-particle system, the 
system normal stresses (pressures) can be obtained from the diagonal spatial components 
of the average energy-momentum tensor (r"), the system energy, from its zeroth diagonal 
component (t°°), and the shear stresses (which are absent for the case of a uniform magnetic 
field) from the off-diagonal spatial components (r*-'). Then, to find the energy density, 
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and pressures of the dense magnetized system we need to calculate the quantum-statistical 
averages of the corresponding components of the energy-momentum tensor of the fermion 
system in the presence of a magnetic field. 

These calculations were carried out long time ago in Ref. j27|, where it was used a QFT 
second-quantized approach. There, a quantum-mechanical average of the energy-momentum 
tensor in the eigen-states of the Dirac equation in the presence of the uniform magnetic field 
was first performed to get the corresponding quantum operator in the occupation-number 
space. The macroscopic stress-energy tensor was then found by averaging its quantum 
operator, in the statistical ensemble using the many-particle density matrix. In this Section 
we perform similar calculations, but using a functional-method approach that makes easier 
to recognize the thermodynamical quantities entering in the final results. Our procedure is 
also different from that of Ref. [^^I in the sense that we will not assume that the fermion 
fields entering in the definitions of the energy and pressures satisfy the classical equation of 
motions (i.e. the Dirac equations for ip and but the functional integrals integrate in all 
field configurations. Then, we keep the terms depend.ng on the Lagrangian derrsity in 
Eqs. dSHD-dSl]), while in Ref. [271 the condition C^= was considered. 

The quantum-statistical average of the energy-momentum tensor is given by 



Tr 



{r'') = ^ (40) 

where 

= [drjd'xKl+r^;^] (41) 
and Z is the partition function of the grand canonical ensemble given by 

Z = Tre-^(^-^^) (42) 

with H denoting the system Hamilonian, N the particle number, /3 the inverse absolute 
temperature, and the chemical potential. The partition function ( H2|) can be written as a 



functional integral 



28| 



Z = j [c/0]eXf / d:'^'^i-'^) (43) 

where £(r, x) is the many-particle system Lagrangian density (i.e. with yU 7^ 0). The 
chemical potential enters in £(r, x) as a shift in the zero component of the electromagnetic 
potential Aq Aq - ^ |29|. 
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A. Energy Density 



To calculate the system energy, (r°j + r°°), we should notice that the external applied 
magnetic field behaves as a classical field in this formalism. Thus, {tIj) = (3V^, being the 
pure Maxwell contribution, with V denoting the system volume. 

To get the fermion contribution to the energy we will make use of the functional integral. 
We start by calculating the quantum-statistical average 



where 



and 



r 



00 



/3 



dr I d^xr'^^iT.x 



(44) 



(45) 



C^ = ^jj [i-f\d° - i/j,) - i-f^d^ - 27^(52 + ieHx^) - i-f^d^ -m]^ (46) 

is the many-particle Lagrangian density. Notice that in (jH]) we are not integrating in the 
photon field A^, since we are taking the Maxwell field as an external classical field. Then, 
to find the fermion contribution to the energy density we have to calculate jy{t^)i with 
given from f lS7|) as 



foo = dr j d^x (#7°aV) + 

dr j d^x'^ [i7°(9° - i7^9^ - ^7^(9^ + ieHx^) - i-f^d^ - m] (47) 
Doing the variable change, r — /3r, in the many-particle partition function, we find 

Z= /'wV][#]e^^o'^-/'^'-^;. (48) 



where 



7°9° + /i7° - ij^d^ - ij^id'^ + ieHx") - ir^d''^ - m 



3q3 



Then 



= I [d^lJ][dii] I (5 dr j d^xC'^ - fi J^' dr J d'x^]'-^^ [ e^I^'^^I^'^ 



-i7°a° 



f3 



(49) 



(50) 



Reversing the variable change (i.e. making /3r — t- r) we obtain 



{dZ/dp) 



(51) 
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where we introduced the grand canonical potential $ = —^InZ, and took into account that 
(A^) = —{d^/dfi)Ty, with N = J d'^xtpj^ip being the particle-number operator. 

Taking into account that the grand canonical potential $ is related to the thermodynamic 
potential by $ = VQ, and that the system entropy is defined as = — (fp)^ ? and the 
particle density by A/" = — we can add the pure Maxwell energy density to Eq. ( ISTll 



1 7t2 

^ '-00 I ~00\ " ~ " . r -n 



e = — {tI' + f""^) = nf + TS + ^,U+— (52) 



to get the system energy density e given by 

^/~00 I ^00\ _ 
B. Longitudinal Pressure 

As in the energy case, to calculate the pressures we will make use of the functional 
integral. For the parallel pressure we start by calculating the quantum-statistical average 

r/3 



(f")-(g + fg) = ^'"""-^'' 7 (53) 



where 

r/3 



^33 



dT j (fxT^^{T,x) (54) 



and is given in (l46!l . 

For the matter field we have specifically 

r/3 



^-33^ ^ immf'M'^i''^'^ ^^^^ 



with given from (138|) as 



+ drj d^x^/j [i7°9° + /i7° - t-f^d^ - i-f\d'^ + ieHx^) - t-f^d^ - m] V (56) 
Now, making the variable change x^ — )■ Lx^ in the partition function, we have 

— jyrjd^.'c'^ (57) 



where 



i7°9° + /i7° - i-f^d^ - n\d^ + ieHx^) - Wd"^ - m 

1j 



^ (58) 



with L a scale factor in the x^, direction and / d^x' = L J d^x ^\ dxs = L J d^x. 
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Then 



L^ = J[d^][(hp]i^LjyT j SxC'^ + L f^dr Jd^x-^'-^^yio'^-I'''^'^'^ (59) 

Reversing the variable change (i.e. making Lx^ — ?■ x^) we obtain 

,00, , (dZ/dL) L (d^\ 



which can be expressed in terms of the thermodynamic potential as 



V 



(r^^) = --Qf (61) 

where V = LA± was considered. 

Following a similar procedure for the pure magnetic contribution and taking into account 
f l28|) . we find, adding the matter and field contributions, 

+ -11) = -^^/-^^ (62) 
Hence, the parallel pressure is given by 

pii = ^n'+-M) = -^f-^ (63) 

C. Transverse Pressure 

To find the fermion contribution to the transverse pressure we start from 

J[dndW:^eIo'^f''^'-^^^^^ 
) = (64) 

where 

rj = dT j d^'xri; (r, x) , 2 = 1,2 (65) 

The explicit form of the transverse diagonal components of the energy-momentum tensor 
f p9|) in the Landau gauge are given by 

f^i = V^27^9V + C^, (66) 

= ?(^7'5' - eH^^x')^/j + C^. (67) 

Apparently, r^^ and r^^ are different, but this is a consequence of the asymmetric Landau 
gauge we are using. Because the magnetic field is along the 0:3 axis, there is an 0(2) 
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symmetry in the xi — X2 plane. Hence, the macroscopic pressures, which are obtained after 
the quantum-statistical average is taken, have to be the same along the xi and X2 directions 
((r^^) = Thus, we can define the transverse pressure as 

1 

2' 



{fi^} = Mfl'} + {rf)) (68) 



with 



r^^ = \j^ dT j (fx^[i-i^d^ + i-i^d^ - eH-i^x^]iP + (69) 



where 

= dr j £xC^ (70) 

Taking into account that in a uniform magnetic field the transverse motion of charged 
fermions is quantized in Landau orbits with radii given in units of the magnetic length 
Iff = 1/ y/eH, we make the variable change 

xj- IhX^, X^ = {Xi,X2) (71) 

in the partition function 

Z = j [<kl)][d'4)]e^'-^ (72) 

where 

S'^ = l]j dr j d=^x^[z7ll ■ d\\ - ilH^-i^d^ - iln^-f^d^ - In^-f^Xi - m]tlj (73) 

Taking the derivative of the partition function Z' with respect to the magnetic length, 
and then reversing the variable change {lu^i 2;^), we obtain 

Ih dZ 



2 dl 



ri^e'^ = (ri^) ■ Z (74) 
From where we have 



Taking into account that 



d _ d (dH\ d _ d d 

+ - 7^ - ^^^H 7T77 (76) 



din din {OIh J dH din OH 
and that V = LA± = Lirljj, we can rewrite ( 1751) as 



-pvn, + pv[H^j (77) 
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Similarly to the longitudinal pressure case, the pure magnetic contribution is 

{rii^) = VP^, 1 = 1,2 (78) 

Finally, we can find the transverse pressure adding the matter contribution ( I77|) and the 
field contribution fl75|) . as 

P± = ^ (^^^^ + r//) = -Qf -HMf + ^ (79) 
where Mf = —{dQf/dH) is the fermion-system magnetization. 

D. Energy-Momentum Tensor Covariant Structure for the Magnetized Fermion 
System 

In analogy to Eq. ( 128|) for the energy-momentum tensor of the magnetic field, we give here 
a covariant decomposition for the energy momentum tensor of the whole system containing 
the matter and field contributions. In this form, we summarize the results for the energy 
density and pressures given in Eqs. fl5^ . (jHS]) and fl7^ . In order to accomplish this goal, we 
define the system thermodynamic potential as the sum of the matter and field contributions 

^ = ^/ + ^ (80) 

Then, we have 

-^{rn = ^V^" + (/^AT + TS)u^u'' + HMr]>'j_\ (81) 

fj V 

where t]'^ was defined in ( 128|1 . and M = —{dVt/dH). 

In the quantum field limit, i.e. when T = fi = H = 0, the only term different from zero 
is the first one in the RHS of ( 18T|) . In that case the system has Lorentz symmetry. If tem- 
perature and/or density are switched on, then the Lorentz symmetry is broken specializing 
a particular reference frame comoving with the medium center of mass and having a four 
velocity = (1, if). This is refiected in the second term of the RHS of ( IHTl) . Finally, when 
there is an external uniform magnetic field acting on the system, the additional symmetry 
breaking 0(3) — )■ 0(2) takes place, and {f^'-') get an anisotropy reflected in the appearance 
of the transverse metric structure rij^ in ( IHTjl . 
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V. EQUATION OF STATE 



As seen from Eqs. (152|) . (l63l) and ( 179!) . the system energy density and pressures depend 
on the fermion system thermodynamic potential Qf. For the dense (/i 7^ 0) fermion system 
under an apphed uniform magnetic field with Lagrangian density (H^ the thermodynamic 
potential is given by (See Appendix [A]) 

— E d{n) dps En + ^In (l + e-^^^"-^)) f 1 + e-^(^"+^)^ 



n=0 



(82) 



where (i(n) = 2— 5„o is the spin degeneracy of the n Landau level, and En = ^Jpl + m? + 2qHn 
is the energy of the particle in the n Landau level with n = 0, 1,2, .... In the bracket of 
the RHS of Eq. fl82l) the first term is the QFT contribution which is independent on the 
temperature T and chemical potential /i, and the second term is the statistical contribution 
depending on these two parameters. 

The particle number density can be obtained from ( 182|) as 



^^ 



my 



2^^^( 



n) 



n=0 



dp3 



I _|_ g/3(en-/i) I _|_ g^(en+M) 



(83) 



We are interested in systems, as neutron stars, where the leading parameter is the fermion 
density /i. Neutron stars cool rapidly through neutrino emission. They can reach in a 
few seconds temperatures T ^ 10^ K fSOi]. Thus, in applications to astrophysical compact 
objects, it is considered that the stellar medium is highly degenerate, so usually neglecting 
the thermal effects. Then, to find the zero temperature limit (/3 — 00) of the thermodynamic 
potential, we will consider only the statistical part in ( l82l) . 



a", 



27r2 ^ ^ ' 



n=0 



dpz^ pI + al - fiJ - a: 



(84) 



where a\ = m? + 2qHn and Uh = -^[(/^^ — m?) /2qH], with /[...] denoting the integer part of 

the argument. The zero temperature limit of the particle number density ( l83l) is 

qH 



n=0 



(85) 



The system magnetization depending on the chemical potential can be found from Eq. 
( ]A22p in the Appendix as 



M 



dH 



J dpsln (1 + 



^^^d{n)J^ dpsln 



-/3(e„-/i) 



I _|_ g-^(en+/i) 
1 



+ 



X _|_ g/3(en-M) 1 -|- g/3(en+M) 



(86) 
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and in the zero-temperature limit it takes the form 

M' = -7r-2T.d{n) r rfp3(^n-/i + ^) (87) 

2tt' ^0 Jo En 

The EoS of the dense magnetized system at zero temperature will be given by the inter- 
relation of the energy density, and the parallel and transverse pressures at zero temperature, 
which can be obtained from (|52|). (|63|). (179|). (|84|). (l85|). and (l87l). as 

e' = Q' + fiAr', Pi\' = -n', p^' = -Q' + H^ (88) 

where is the general thermodynamic potential of the system, given by 0° = fij + 

In the case of quark matter, the asymptotically-free phase of quarks will form a perturba- 
tive vacuum (inside a bag) which is immersed in the nonperturbative vacuum. This scenario 



is what is called the MIT bag model SjJ- The creation of the "bag" costs free energy. Then, 
in the energy density the energy difference between the perturbative vacuum and the 
true one should be added. Essentially, that is the bag constant B characterizing a constant 
energy per unit volume associated to the region where the quarks live. From the point of 
view of the pressure, B can be interpreted as an inward pressure needed to confine the quarks 
into the bag. Hence, in the case of quark matter, Eqs. (|88ll get extra terms depending on 
B, 



e' = n^+^Af^ + ^ + B, p||0 = -f]0 p/ = _f]0 /^M° + ^-i? (89) 

■'2 •'2 ■'2 



Because the relative sign between B and the magnetic energy /2 term is not the same 
for the parallel and transverse pressures, the pure magnetic energy contribution cannot be 
absorbed by the vacuum energy B. Taking into account that the magnetic field varies in 
several orders from the inner core to the star surface, the term applies a tremendous 

extra pressure on the quark matter, but because of the anisotropy between the longitudinal 
and transverse directions with respect to the field alignment, the pressure coming from 
is negative on the parallel pressure, while on the transverse pressure it is positive. 

The phenomenological parameter B is estimated taking into account the underlaying 
dynamics and external conditions (as temperature and density) of the system SsJ, but it 
cannot be calculated from first principles due to our present limitations in dealing with 
non-perturbative QCD. In our analysis of the magnetic field dependence of the equation of 
state, we will only consider the free-fermion case, hence we take i? = 0. 
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(a) (b) 



FIG. 2. (a) System energy density vs magnetic field for /i = 'iOOMeV. (b) System pressures 
(parallel pressure represented in solid line, and transverse pressure in dash line) vs magnetic field 
for fi = AOOMeV. 



In Fig. 2, we are showing how the energy density and pressures change with the magnetic 
field at fixed /i. We found that the energy density and transverse pressure increase with the 
magnetic field, while the parallel pressure decreases. Also we have that, in our free-fermion 
model, for field strengths close to lO^^G the parallel pressure becomes negative. Hence, field 
strengths of that order can produce strong instabilities in the star's structure. It will be 



interesting to investigate this issue in more realistic models that include partic 
The splitting between the two pressures relative to their weak field value 



e interactions. 



LJ (QQ) 

" I p||0(e5 < /i2) I ^ ^ 

is given in Fig. 3. There, we can see that for /i = AOOMeV, fields H > lO^^G have splitting 
rates A > 10. This result indicates that for the field-strength range that can take place in 
the inner core of magnetars, for example, the pressure anisotropy can play an important role 
in the star structure and geometry. A criterium to separate the isotropic regime {p±q = P\\q) 
from the anisotropic one can be 

A ~ C(l) (91) 

For the density value previously considered (/i = 400MeV), we have a threshold field of 
order lO^^G. 



Notice in Fig. 2 that at weak field {eH ^ /i^) the two pressures coincide. 
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FIG. 3. Splitting coefficient A vs magnetic field for fx = 4Q0MeV 




log^„( qH / iL^' ) log^„( qH / m? ) log^„( qH / jl^ ) 



(a) (b) (c) 



FIG. 4. (a) Ratios of system/matter energy density, (b) parallel pressure, and (c) transverse 
pressure vs magnetic field for fi = AOOMeV. 

To determine for what field values the pure magnetic contribution to the energy density 
and pressures becomes important, we plotted the ratio between the total energy (i.e. the 
one having the matter and field contributions) and the matter energy density in Fig. 4- (a); 
that of the total parallel pressure to the parallel matter pressure in Fig. 4-(b); and the 
corresponding one to the transverse pressure in Fig. 4-(c). From those graphs we have 
that the pure magnetic contribution becomes significant for field strengths between one and 
two orders smaller than //^. Thus, for n — AOOMeV, fields larger than lO^^G will make a 
significant contribution in the system energy and pressures through the term. Hence, 

24 




FIG. 5. System magnetization M times the magnetic field H vs the magnetic field for ji = 400Mel^. 

in astrophysical applications where the field strength in the inner core of compact objects 
can be large enough, the pure field contribution should be consider on equal footing as the 
matter contribution in determining the parameters of the EoS. It is worth noticing that the 
pure magnetic contribution becomes significant for fields of the order of the threshold field 
for the isotropic-anisotropic transition. This is not a coincidence, but a consequence of the 
fact that the main contribution to the pressure splitting 

= ifM° + i7^ (92) 

comes from their magnetic term H^. The term depending on the magnetization M° has an 
oscillatory behavior due to the Haas-van Alphen oscillations appearing at low temperature 
in degenerate fermionic systems [ssl (See Fig. 5). Nevertheless, the Haas- van Alphen 
oscillations are not noticeable in Fig. 3 because their amplitudes are much smaller than H, 
so they can be neglected in fl92l) . 

The threshold value of ~ lO^^G we found, is model dependent. It is applicable to systems 
of free fermions under certain parameter values. For other system, such as cold dense quark 
systems with superconducting gaps, the corresponding threshold field should be determined. 

Finally, in Fig. 6 the system EoS is plotted . There, the variation of the energy density 
vs the parallel pressure is given in Fig. 6-(a), and with respect to the transverse pressure 
in Fig. 6-(b), for a fixed density (/x = 400MeV^) and a variable magnetic field. Due to the 
pressure anisotropy, the EoS in this case should be given by a curve in a three dimensional 
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FIG. 6. (a) System energy density minus the bag constant B vs the parahel pressure and (b) vs 
the transverse pressure, for magnetic fields in the range < qH/ji^ < 1 and /x = AOOMeV. 

representation with axes {e,p\\,p±). In Fig. 6 we give the projections of that curve in the 
two planes, both including the e axis. 




VI. CONCLUDING REMARKS 

In this paper we have shown that a relativistic dense fermion gas under a sufficiently 
strong magnetic field has a highly anisotropic EoS. For field strengths about two orders 
smaller than /x^, the anisotropy effects begin to emerge, and once the field reaches values 
one order smaller than /x^, they cannot be neglected anymore, since the splitting is ten times 
the value of the pressure at zero field. The splitting of the pressure in two terms, one along 
the field direction (the parallel pressure) and the other, perpendicular to the field direction 
(the transverse pressure) , should be taken into account in astrophysical considerations when 
studying compact objects that exhibit strong magnetic fields, as it may affect the structure 
and geometry of the star. 

At strong magnetic fields {H ~ /x^/lO), the pure magnetic energy contribution £m, as well 
as the magnetic pressures p\\j^ and p±M, were found to be as important as the corresponding 
matter contributions Ef, py^ and p±f (see Fig. 5 for details). Therefore, a sufficiently high 
magnetic field may actually infiuence the EoS in two different ways: by modifying the 
matter contributions to the energy density and pressures, and as importantly, through the 
pure magnetic (coming from the Maxwell term in the original Lagrangian) contribution to 
the energy and pressures. The order of the fields required for the two contributions to 
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become comparable was found to be the same as the field strength needed for the pressure 
anisotropy to be relevant. We call attention that the magnetic contribution into the EoS has 
)een proved to have significant effects in the magnetohydro dynamics of quark-gluon plasma 



34| 



Our results are valid for relativistic systems of fermions in the presence of a uniform 
and constant magnetic field. These kinds of systems have been considered in the physics 
of neutron stars, as well as quark stars. At present, it is known that if the density is high 
enough to produce a quark deconfined phase in the star core, the state that minimizes the 



energy of such system will be a color superconductor [35|. There exist already several papers 
where the effect of the color superconducting gap was considered in the energy density and 
pressure of highly dense quark matter [36|] . 

In spin- zero color superconductivity, the color condensates in the CFL, as well as in the 
2SC phases, have non-zero electric charge. Then, we could expect that the photon acquires 
a Meissner mass which produces the screening of a weak magnetic field (the well known 
phenomenon of Meissner effect). Nevertheless, in the spin-zero color superconductor the 
conventional electromagnetic field is not an eigenfield, but it is mixed with the S'^'-gluon 
to form a long-range field that becomes the in-medium electromagnetic field (i.e. 



the so called "rotated" electromagnetic field) 



37| . In a series of papers 
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there has been 



shown that the color-superconducting properties of a three-fiavor massless quark system are 
substantially affected by the penetrating "rotated" magnetic field and as a consequence, 
a new phase, called Magnetic Color Flavor Locked (MCFL) phase, takes place. In the 
MCFL phase the gaps that receive contributions from pairs of charged quarks get reinforced 
at very strong fields producing a sizable splitting as compared with the gaps that only 
get contribution from pairs of neutral quarks. As the field decreases, the gaps become 



oscillating functions of the magnetic field 39(|, a phenomenon associated with the known 



Haas- van Alphen oscillations appearing in magnetized systems [33 1. 

Therefore, a realistic study of the EoS of stellar magnetized dense quark matter should 
be formalized in the color superconducting MCFL phase. A preliminary study already 
exists, where the magnetic field effect was considered only in the quark particle spectrum, 
but no in the gap, which was taken as a fixed parameter. In a future publication 4l|], a 



self-consistent investigation of the EoS in the MCFL phase, taking into account the gap 
equations depending on the magnetic field, will be presented. 
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Appendix A: Thermodynamic Potential for a Dense-Magnetized Fermion System 



The thermodynamic potential, Qf, in the presence of a constant and uniform magnetic 
field has been previously calculated using different methods at finite temperature and/or 



chemical potential. Fo r example, using the Schwinger proper time method 
calculated in Refs. 



421 ]. Qf was 



46|, and 



431 and Ml, at T ^ 0, and introducing a chemical potentia 

n 

respectively. In the Furry picture [45'], Qf was calculated at T 7^ and 7^ in 
using the worldline method at T 7^ in 47]. 

Here we present in detail the calculation of the thermodynamic potential in the presence 
of a constant and uniform magnetic field at T 7^ and /i 7^ usin g R itus's method. 
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and later on 



This approach was originally developed for charged spin-1/2 particles 
extended to charged spin-1 bosons 49|. Recently, it has been implemented for the case of 
spin-1/2 in an inhomogeneous magnetic field in reduced dimensions jsO]. This approach 
is based on a Fourier-like transformation performed by the eigenfunction matrices Ep{x) 
which are associated to the wave functions of the asymptotic states of charged fermions in 
a uniform magnetic field. The Ep{x) functions play the role in the magnetized medium of 
the usual plane-wave (Fourier) functions e*^^ at zero field. The advantage of this method 
is that the field-dependent fermion Green function is diagonalized in momentum space, so 
having a similar form to that in the free space. Hence, this formalism is very convenient to 
implement the statistical sum by the imaginary time procedure needed to describe systems 
at finite temperature and density . Also, the obtained Green function in momentum space 
explicitly depends on the Landau levels. This last result makes it particularly convenient to 
be used in the strong-field approximation, where one can constraint the calculations to the 
contribution of the lowest Landau level (LLL) in the particle spectrum. 

Ritus's method has been successfully used in the context of chiral symmetry breaking in 
a magnetic field [si], as well as in magnetized color superconductivity 38], 52]. 

The thermodynamic potential of the magnetized dense system at finite temperature 
Qf{H,fi,T) is given by 



nf{H,fx,T) 



V 



(Al) 
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where ^{H,fi,T) is the grand canonical potential (in functional terminology, the effective 
action in the presence of an external magnetic field at finite temperature and density), which 
is given in terms of the inverse fermion propagator as 

<^>{H,fi,T) = ^TrlnZ = ^Tr In G-\x, x'), (A2) 

with the trace and logarithm taken in a functional sense, and G^^{x,x') being the fermion 
inverse propagator in space representation. To make the transformation to momentum space, 
because of the dependence of G~^{x, x') on the electromagnetic potential of the external field 
j^ext _ }jxiSfj_2, it is convenient to use the Ritus's transformation 

G-\x,x') = f^j^,K{x)mGT\v)E[{x') (A3) 

where = "^S/^o ^^f^p^! and we introduced the Ritus' transformation functions 

E^ix) = E;{x)^{+) + E;{x)^{-) (A4) 

with 

A(±) = (A5) 

representing the spin up (+) and down (— ) projectors, and Ep~{x) are the corresponding 
eigenfunctions 

E+{x) = iVie-^(Po^°+f^^'+f^^')A(p), 
Ep{x) = A^,_ie~*(P«^°+P^^'+f^^')A-i(p) (A6) 

with normalization constant A*"; = (AneHy^'^/y/ll, Di(p) denoting the parabolic cylinder 
functions of argument p = \/2eH[xi —p2/eH), and index given by the Landau level numbers 
/ = 0,1,2,.... 



The Ep functions satisfy the orthogonality condition 53 1 



/ d'xE'^{x)E^;,{x) = (27r)^5W(p - p'Ml) , (A7) 



with Ep = -f%El,y-f°, U{1) = A(+)5'o + /(I - S^^), and S^^\p - p') = (5"'5(po - Po)SiP2 - 

The spin structure of the Ep functions is essential to satisfy the eigenvalue equations 

{U-^)Elix) = El{x)i^-p) , (A8) 
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with p^' = 0, -V2eHl,p^). 

Using (jMl) and (jA8l), the inverse propagator in momentum representation Gi (p) ap- 



pearing in (IA3p . is found from 



= (2vr)^5(^)(p-p')n(/)Gr^(p) 



where 



and 



n^, = id,, - eA^ 



Gi (p) = [p* ■ 7 - m] 



(A9) 



(AlO) 



(All) 



with pI = (ip'^ - fi, 0, V2eHl,p^) 

Substituting (lASp in (1A2|) . taking the functional trace, and using he orthogonahty con- 
dition ( ]A7p . we obtain 

74, 



$(ff,/i,T) = -trln / dx I dx'6\x,x')f^^,El{x)U{l)Gr\p)Ep{^') 



% 

T 

?3/ 



trln^U{l)GY\p)d' 



p 



(A12) 



In f lA12p . d^p = dpodp2dps, and tr denotes the remaining spinorial trace. Now, taking into 
account that 



^li^) = TTT^ / dx4 / dX2dX3 

^ [ZTC)-^ jo J-oo 



(A13) 



and because {p) does not depend on p2, for the integration in p2 in ( ]A12I) we have 



°° dp2 f°° dp2 _iP2n. 1 - eH 



-e 



lpi=0 



dxi 



l-oo 271 J-oo 27r 

Substituting (lAlSp and (]A14|) in (]A12I) . we obtain in Euchdean space {po — )• 2^4) 



(A14) 



00 „ 

<l>{H,fi,T) = -eHVptrlnY, / 



dpdpz ^n\r<-i(-\ 



(A15) 



Taking into account that n(/) separates the / = Landau level from the rest, substi- 
tuting flAlSP in flAll) . and because of the identity trlnO = IndetO, we obtain for the 
thermodynamic potential 



nf{H,fi,T) = -eH 



^^^^^^lndetGo-(p) + 2f: ^^^^^^ 



00 (27r 



1=1 



{2nf 



In det G," 



(A16) 
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The fermion system at finite temperature can be described by taking the discretization 
of the fourth momentum following Matsubara's procedure 



dp4 1 ^ 



P. = ^^^, n = 0,±l,±2, 



(A17) 



Taking into account (1A11I) . after taking the trace in (]A16I) we obtain 

dps 



QfiH, /z, T) = r ^ E(2 - '^^o) E In [(P4 + + e 



where 



ef =pI + 2\eH\l + 



(A18) 



(A19) 



Notice that in this approach the sum in the p^ term, which is obtained in ( lAlSp . is formally 
similar to that appearing in the free-particle thermodynamic potential (i.e. at if = and 
/i 7^ 0) 28|. After summing in p^ we get 



nf{H,f^,T) 



^ ^eHdil) 
^^P3E^^ 



(A20) 



The ratio with d{l) = 2 — is the density of states per Landau level. The factor 

d{l) is the spin degeneracy of Landau levels with / 7^ 0. 

The thermodynamic potential flA20p has two contributions. One that does not depend 
on the temperature and chemical potential Q'^^^, and the statistical one given 
respectively by 



nf^-iH, T, /.) = J_l dp, E dH) In (1 + e-^^-'^)) (l + 



— 2 / dps^d{l)ei 



47r2 



1=0 



(A22) 



As known, Vfi^^jH) has non- field-dependent ultraviolet divergencies that should be renor- 



malized (see Ref. 



54| for a detailed renormalization procedure of this term). After renor- 



malization the well-known Schwinger expression 



42| 



q9/^{H) = --^ r ^expi-m's) ( esH coth{esH) - 1 - ]:{esH)A (A23) 
■' Svr^ Jo s-^ \ 3 / 

is found. In the calculations of the energy density and pressures done in Sec. V we considered 
only the Q^'^^'^ contribution, since for astrophysical applications one always has /x^ ^ eH, 
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and the leading contribution oiQf will come from Q 
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